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ABSTRACT 

The proton burning process p + p — > d + e + + v e , important for the 
stellar evolution of main-sequence stars of mass equal to or less than that 
of the Sun, is computed in effective field theory using chiral perturbation 
expansion to the next-to-next-to-leading chiral order. This represents a 
model-independent calculation consistent with low-energy effective theory 
of QCD comparable in accuracy to the radiative np capture at thermal 
energy previously calculated by first using very accurate two-nucleon 
wavefunctions backed up by an effective field theory technique with a 
finite cut-off. The result obtained thereby is found to support within 
theoretical uncertainties the previous calculation of the same process by 
Bahcall and his co-workers. 
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1. Introduction 

The proton fusion reaction 

p + p — > d + e + + z/ e (1) 

which plays an important role for stellar evolution and - as the dominant neutrino 
source - for the solar-neutrino problem, has quite a long history of investigation. 
Indeed the reaction rate of this process (hereafter called the pp rate) was first 
calculated by Bethe and Critchfield ( |1938| ). Salpeter ( |1952| ) recalculated the pp rate 
using the effective range approximation and argued that the relevant nuclear matrix 
element squared could be estimated with an accuracy of the ~5 % level. (The pp 
rate itself was subject to much larger uncertainty, ~20 %, because of the limited 
precision with which the Fermi coupling constant was known at that time.) Bahcall 
and May ( |1969|) examined the dependence of the pp rate on explicit forms of the 



two-nucleon wavefunctions generated by two-parameter nuclear potentials of various 
forms adjusted so as to reproduce the scattering length and effective range (for the 
pp channel) and the low-energy properties of the deuteron (for the np channel). The 
pp rate was found to vary by ~ 1.5 % corresponding to the changes in the deuteron 
wavefunction, and by ~ 1.2 % due to the change in the pp wavefunction. The 
most updated work along this line was done by Kamionkowski and Bahcall (|1994|) 
employing deuteron wavefunctions obtained from much more accurate potentials 
such as the Argonne vu, t>i 8 , Urbana vu, super-soft-core (SSC) and Reid soft-core 
potentials. Changes in the pp rates arising from the different potentials were found 
to be ~ 1 %. Thus it seems that the presently available calculated pp rate is 
robust and needs no further scrutiny, the famous solar neutrino problem remaining 
unresolved from this angle and hence persisting as one of the outstanding unsolved 
problems in astrophysics (Bahcall |1997| ). 

There are however two reasons for revisiting this issue. One is that while 
the calculated pp rate seems to have converged to a "canonical" value given in 
(Kamionkowski & Bahcall |1994|, hereafter KB), there lingers the unsettling feeling 



that the strong interaction involved in nuclear physics of the two-nucleon systems 
is infested with uncontrollable uncertainties associated with model dependence in 
the treatment, making it difficult to assess the accuracy achieved. Thus it is not 
unexpected that this canonical value will be - as has been in the past - challenged. 
Indeed it has recently been argued by Ivanov et al. (|1997|) that the nonrelativistic 
potential models used in the previous works could be seriously in error. They show 
that in their version of a relativistic field theory model, the pp rate comes out to 
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be as big as 2.9 times the previous estimates.^ 5 ! Should their new result turn out 



to be correct, it would have profound consequences on theories of stellar evolution 
in general and on the solar neutrino problem in particular. In a nutshell, the 
issue comes down to whether or not a more general framework such as relativistic 
field theory would invalidate the calculation made in the traditional nonrelativistic 
potential models. The claim of the authors in (Ivanov et al. |1997|) is that it indeed 



does. Our aim is to address this issue using a low-energy effective field theory of 
QCD that has found a quantitative success in other nuclear processes. 

The second reason is really more theoretical, independent of the above 
important astrophysical issue. Along with the thermal np capture, the proton 
fusion process is the simplest nuclear process amenable to an accurate calculation 
- something rare in hadronic physics - and it is of interest on its own to test how 
well a calculation faithful to a "first-principle approach" can tackle this problem. 
In particular, we are interested in checking how accurately the effective field theory 
approach, found to be stunningly successful for the np capture n + p — > d + 7, 
low-energy NN scattering and static properties of the deuteron (Park, Min, & Rho 
1995| ; Park et al. |1997b|) , fares with the proton fusion problem, the weak interaction 



sector of the Standard Model. 

The strategy we shall adopt here is quite close to that used for the np 
capture process (Park et al. |1995| ). We shall use chiral perturbation theory to the 
next-to-next-to-leading order (NNLO) in chiral counting; as defined precisely below, 
this corresponds to 0(Q 3 ) relative to the leading-order term. This is roughly the 
same order as considered for the np capture. However the relative importance of 
terms of various chiral orders is somewhat different here. As we explain later, in 
the present case, the corrections to the leading order are not suppressed by what is 
called the "chiral filter" W*\ and so the accuracy with which these can be calculated 



# 6 They use a procedure that seems to disagree with other physical properties of low-energy pp 
systems, as was pointed out by Bahcall and Kamionkowski (1997). It has also been pointed out 
by DegPInnocenti et al. ( 199S|) that such a large deviation from the value used by KB would be 



inconsistent with hclioscismology in the Sun. 



# 7 The chiral filter phenomenon is explained in detail in (Park, Min, & Rho, 1993). Crudely stated, 
it refers to the general feature that whenever one soft-pion exchange is allowed by kinematics and 
selection rules, it should give a dominant contribution with higher-order (or shorter-range) terms 
strongly suppressed. A corollary to this is that whenever one soft-pion exchange is not allowed, 
all higher-order terms can be important, making chiral perturbation calculation generically less 
powerful. 
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is not as good as in the np capture case. Even so, using the argument developed 
in (Park et al. |1997b|) , we shall suggest that the procedure used here provides a 
model-independent result in the same sense as in (Park et al. |1995| , Park et al. 
I997T5D . 



To streamline the presentation, we first give our result and then discuss 
(as briefly as possible) how we arrive at it in the rest of the paper. Apart from 
the meson-exchange contributions which are of order of 0(Q 3 ) and which for the 
reason mentioned above and further stressed in our concluding section, are the main 
uncertainty to the order considered, our chiral perturbation theory result in terms of 
the reduced matrix element A defined in (Bahcall & May |1969[ ) is 

A 2 xPT = (1 ± 0.003) x 6.93 (2) 

where the uncertainty is due to experimental errors]^ The above result is to be 
compared with the value obtained by Kamionkowski and Bahcall (KB) 

A| B = (l ± 0.0021^) x 6.92 . (3) 

As we shall explain, there are some differences in details between our calculational 
framework and that of KB, but our final numerical result is in good agreement with 
that of KB and disagrees with that of Ivanov et al. (|1997|) . 



The paper is organized as follows. In Section |2|, our strategy for carrying out 
a chiral perturbation calculation for two-nucleon systems is outlined. Our approach 
here is similar to the one used in the previous calculation of the np capture process. 
We shall sketch a justification of this approach from the standpoint of low-energy 
effective field theory of QCD (with concrete supporting evidence summarized in 
Subsection ^2|) . Section || describes chiral counting of the terms appearing in the 
relevant weak current. In Section f|, the wavefunctions for the initial pp state and 
the final d state are specified. Our numerical results are given in Section [5], and a 
brief discussion including a comment on the main uncertainty in the calculation is 
given in Section H 



2. Effective Field Theory for Nuclei 

Before presenting our calculation, we sketch the chain of reasoning that leads 
us to believe that the result we obtained is model- independent. The argument is 



#s Our theoretical uncertainty is, if very conservatively estimated, about 0.1 %. 
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essentially identical to that used for thermal np capture, and therefore we review 



briefly the key ingredients of the calculations in (Park et al. [I995| ; Park et al. |1997b|) 
for the process 

n + p — > d + 7. (4) 
The basic premise of our approach is an application of the "Weinberg theorem" 



(see, e.g., Weinberg |1996|) , which for the case in hand implies that at very low energies 
at which the process @ is probed, QCD can be described by chiral perturbation 
theory. In setting up chiral expansion for two-nucleon systems we are interested in, 
we follow Weinberg ( |1990|) and separate the relevant transition amplitude into the 
"reducible" and "irreducible" terms. Bound states and resonances involve reducible 
graphs that are kinematically enhanced. To incorporate this feature, we need to 
sum up a given set of reducible terms to all orders, a procedure equivalent to solving 
the Schrodinger or Lippmann-Schwinger equation. The irreducible graphs that 
appear as kernels in the latter are calculated to a specified order in chiral expansion. 
How this can be done in a systematic way consistent with chiral counting has been 
extensively discussed in the literature (Ordonez, Ray, & van Kolck |1994j ; Kaplan, 



Savage, & Wise gggg ; Luke & Manohar |1997| ; Beane, Cohen, & Phillips g997D . In 



(Park et al. |1995|) , we have written the relevant Ml matrix element for the process 
(H) with the EM current given by the irreducible graphs calculated to NNLO and 
accurate wavefunctions for the initial and final states, namely the Argonne vig wave 



function. The Argonne vi$ potential (Wiringa, Stoks, & Schiavilla |1995| ) provides an 



impressive description of all two-nucleon data at low energies up to ~ 350 MeV. So 
the wavefunctions obtained with it correspond to summing both the irreducible and 
reducible graphs to all orders. It seems then that the procedure used in (Park et 



al. |1995|) is not quite consistent since the current is computed only to NNLO. We 
argued however in (Park et al. |1995|) that this is consistent with the chiral counting 
to the chiral order considered and the error made in the terms of higher order than 
that taken in the current cannot be important. We have subsequently demonstrated 
(Park et al. |1997b|) the correctness of this conjecture by showing in a cut-off field 



theory that a systematic summation of the reducible diagrams with the irreducible 
graphs appearing as a kernel - which accurately describes all static properties of 
the deuteron and NN scattering lengths - gives precisely the same result for the 
leading chiral order Ml matrix element as that obtained with the Argonne v is wave 
function. Now, as we know, the ratios of subleading chiral terms to the leading term 
are almost model-independent and insensitive to the renormalization schemes used, 
so these ratios can be calculated with reasonable accuracy. Thus the result obtained 
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in (Park et al. [T555D , 

< rX pT = 334±3 mb ^ 

(where the error represents our ignorance regarding the short-distance part of the 
interaction in chiral perturbation theory) which agrees well with the experiment, 

a ex P = 334.2 ± 0.5 mb, (6) 

can be taken as a first-principle calculation as argued in (Park et al. |1997b|) . 



For the weak process ([!]) we will follow the same procedure and, to supply 
support for it, carry out an effective field theory calculation up to NNLO with the 
use of the same cut-off scheme as in (Park et al. |1997b| ), where the leading-order Ml 
matrix element was evaluated. 



3. Chiral Counting and the Weak Current 



As in KB (|1994 ), we shall include the vacuum polarization effect pertinent to 



the initial pp channel as well as the effect of meson-exchange currents. In addition, 
we include the effect of the two-photon exchange Coulomb interaction for the pp 
channel. 



3.1. Chiral counting 



The core temperature of the Sun is believed to be about T c = 15.5 x 10 6 K, 



for which the corresponding nucleon momentum is p ~ y2m p Tc ~ 1.6 MeV. We 
are therefore interested in the process ([!]) at very low energies, a situation similar 
to the thermal np capture process. This is an arena where the low-energy effective 
field theory can be powerful. At these low energies, the reaction is dominated 
by the 1 Sq — > d transition, involving only the transition operators that carry 
orbital angular momentum L = (0, 2), spin S = 1 and total angular momentum 
J — 1. The momentum carried by the leptons, k = p e +p u , is also very small, 



< Ifcl < 



m e are the mass of proton, deuteron and electron, respectively. We can therefore 
safely ignore terms of (9(|/;;| 2 ) and work only to These selection rules will be 

implicit throughout our subsequent arguments 



i 

2 



0.78 MeV, where m p , rrid and 
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To carry out the procedure explained in Section 0, we need to specify a counting 
rule for irreducible terms for the electroweak hadronic current J M ' _ = J^' 1 — zJ 7 '' 2 , 
where 



This rule has been extensively discussed in (Park et al. |1995| ; Park et al. |1993| ; Park, 



Jung, & Min |1997a]) , so we shall simply summarize it here. 

Let Q denote the typical momentum scale probed in the process and A x be the 
chiral scale ~ mjv- The chiral expansion is made in powers of Q/A x . A Feynman 
diagram for an A-nucleon process (Weinberg |199(J| ) goes as Q u . Now for a Feynman 
graph with L loops, N E external fields and C separately connected pieces, the 
exponent v is given by 

v = 4 - A - 2C + 2L - N E + ^ v% , with iy l = d i + e i + — - 2, (8) 

i " 

where rii is the number of nucleon lines, di the number of derivatives or powers of 
m n and the number of external fields at the i-th vertex. Since z/j is defined so that 
chiral symmetry guarantees > (Weinberg |1990| ) even in the presence of external 



fields (Rho |199"l|) , the effective Lagrangian can be ordered according to uf 

C eff = C + C 1 + C 2 + ---, (9) 
where C n has vertices with z/j = n. 



3.2. Current 



Let us now apply this counting rule to our case. The leading contribution 
arises from a one-body diagram with v = uq = —3 (i.e., L = 0, C = A = 2, N E = 1 
and Vi = 0), and the next-to- leading order [y = u + 1) also comes from a one-body 
diagram involving a ^ = 1 vertex. These one-body operators are given by 

W-£*H,!±? + &fx,). do, 




with pi (p() the incoming (outgoing) momentum of the Z-th nucleon, — 1.2601(25) 
(Particle Data Group 1996p , \iy = \i p — \i n ~ 4.7059 the nucleon isovector 
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magnetic moment and the nucleon mass. Of these one-body operators only the 
Gamow- Teller term contributes at the leading order, while the weak magnetism 
(WM) term that comes at the next order is further suppressed kinematically for the 
small lepton momentum \k\ <C Q. The others do not contribute by the selection rule. 
We will keep the WM term but neglect all other higher chiral order fc-dependent 
terms. Thus the one-body current we calculate is 

J\b = ^gt + -Avm + ' ' ' > (12) 

with 

^gt = -9aY, t i~^i, (13) 

I 

4m = T^-kx^nai. (14) 

To proceed to higher orders, it is convenient to do chiral counting relative to the 
leading order. Expressed in this way, the leading order one-body current eq.(|12D is 
taken to be ~ 0(1), and the next-to-leading order term is ~ 0(Q). The two-body 
exchange current can enter at the next-to-next-to- leading order (NNLO), that is 
at 0(Q 2 ). Because of the selection rule, however, there is no term of this order 
contributing to the process ([I]). Hence the two-body exchange current at tree order 
enters at 0(Q 3 ). Note the difference in the roles of the exchange currents between 
the present case and the np capture process. In the latter, the exchange-current 
(2-body) Ml operator is ~ 0(Q) relative to the leading one-body Ml operator. This 
is due to the fact that the one-body Ml operator is formally suppressed by one order 
of chiral counting, so that the one-soft-pion exchange two-body current is enhanced 
relative to the one-body operator. This is just the chiral filter enhancement referred 
to above. In the Gamow- Teller case, the chiral filter mechanism is not operative, 
making the leading two-body correction come at 0(Q 3 ). This renders evaluation of 
the two-body terms less precise although their overall importance relative to the 
one-body term should be considered diminished here. We will return to this matter 
later. 

Calculation of the meson exchange Gamow- Teller operator is discussed in 
(Park et al. |1997a|) . Here we skip the details and simply quote the result. The 
leading two-body contributions of 0(Q 3 ) come from the one-pion exchange and 
four- Fermi contact graphs with L = 0, C = A — 1 = 1 and one z/j = 1 interaction 
(y = vq + 3). The relevant part of the Lagrangian with the z/j = 1 interaction is of 
the form (Bernard, Kaiser, & Meissner |1997| ; Cohen et al. |1996|) 



_ / V » V V — Q> tv / 1 \ \ 

& = B( D^D V + 4c 3 iA • iA + (2c A + — j [S* 1 , S»] [iA„ »A„] j B 
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4idi BS -ABBB + 2id 2 e abc e„ uXS v ll A v > a BS X T b B BS s r c B + ■■■ 



(15) 



where £0123 = +1, and the definitions of the "covariant derivatives" _D M , A M and the 
spin operator can be found in (Bernard et al. |1997| ; Cohen et al. |1996| ). There are 
no vector-current contributions and the surviving two-body Gamow- Teller operator 
in momentum space is 



A a 



9A 



2m N f* ml + q 2 { 2 ®' 
+ c 3 qq ■ (r^a e + T*a e ) + (c 4 + ^ r° q x (a® x q) J 



9a 



(16) 



with p = |(pi — p 2 )> Pi = + /tt — 93 MeV is the pion decay constant and 
is the pion mass. We have defined 0©,e,® by 



©,e 



and similarly for e <~ . The dimensionless parameters q and dj are defined as 



1 



■Cj and dj 



9a * 
n d i- 



;i7) 



;is) 



Among the l/m^ corrections, terms proportional to P = p 1 + p 2 = —\k have been 
dropped. The p-dependent terms, which come from ^-^BD 2 B in eq. (|T5|) , shall be 
referred to as the "kinetic term". The low-energy constants c 3 and C4 have been 
determined from pion-nucleon experiments by Bernard et al. ( 1997 ): 



c 3 
c 4 



-5.29 ±0.25 GeV" 1 , 
3.63 ±0.10 GeV" 1 . 



(19) 



They are almost completely saturated by the resonance-exchange contributions 



„Res 



Jtes 



4 + <4 + cf 



-3.83 - 1.40 - 0.06) GeV" 1 = -5.29 GeV" 1 , 



cf + 4 + cf = (1.92 + 1-63 + 0.12) GeV" 1 = 3.67 GeV" 1 



(20) 



where the superscripts A, p, S and R denote the contributions from the exchange of 
A(1232), p-mesons, scalar-mesons and the Roper resonance, respectively. In terms 
of the dimensionless constants c 3 and c 4 , eq. (P~8D gives us 



c 3 ~ -4.97 ±0.23, 
c 4 ~ 3.41 ±0.09. 



(21) 
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There are no low-energy data to determine the constants di and d 2 - However, it 
was argued in (Park et al. |1995| ) that the zero-range operators corresponding to the 
d\ and d 2 terms should be killed by the short-range correlation. We shall justify 
this argument here by showing that for "natural" coefficients in the Lagrangian (as 
required by chiral symmetry) their matrix elements are indeed negligible. 



Finally, when Fourier-transformed, the two-body current in coordinate space 



reads 



J. 



2B> 



+ 
+ 



2m N fl \ 2m\r 



+ c 3 I f tJ y 2 (m n r) + —y^m^r) 



+ ( c 4 + - I r" I -r lJ y 2 (m„r) + 



25 ij 



-yoim^r) a; 



9a 



2rriNf% 



1 



di + -c 3 ) (r^ + t»? q ) + ( -2d 2 + -64 



a -n 



(22) 



where the ellipses denote irrelevant terms and the Yukawa functions, y's, are defined 
as 

e -x 

Airx" 1 



3/o W 



:i+x) 



Attx 
, 3 3 

1 + - + -5 



d 



X 



x- 



Anx 



did 

x- — -i-yo{x). 



(23) 



4. Wavef unctions 

In the scheme we are using, it is essential to have accurate wavefunctions 
for the initial and final states so that the leading single-particle Gamow- Teller 
matrix element is accurately given. In (Park et al. 1995Q , we have shown that the 



single-particle Ml matrix element for the np capture process (§) is accurately given 
by the wavefunctions obtained with the Argonne v±s potential. The single-particle 
Gamow- Teller operator is identical in structure to the isovector Ml operator, so 
provided that electromagnetic effects are suitably taken into account, we expect 
that the Gamow- Teller matrix element for the proton fusion process (|lj) will also 
be given very accurately by the wavefunctions of the Argonne Vx$ potential. Once 
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this argument is justified, then we can conclude immediately that the answer 
was essentially given by Bahcall and collaborators since they used, among others, 
the Argonne vi& potential in calculating the Gamow- Teller matrix element. One 
difference in their calculation from ours is that in their treatment, the initial 
and final wavefunctions were calculated asymmetrically with the initial scattering 
wavefunction treated approximately while the final deuteron state treated "exactly" . 
Our strategy in the chiral perturbation approach requires that they be treated with 
the same "exact" potential. Numerically, however, their approximation for the initial 
state should be quite good, so it is not unreasonable to expect that not much will be 
gained by the symmetric treatment. We confirm this. 

We now turn to discussion of the wavefunctions used in our calculation. 



4.1. The deuteron wavefunction 



The deuteron wavefunction obtained with the Argonne fig potential (referred 
to as the "Argonne potential" in what follows) is well-known, so without going into 
any details, we will merely define the notations for later use. 



The deuteron wavefunction will be written in the conventional form 

A s 



Airr 



Ud{r) H 7 =-w d (r) 



'Si), 



(24) 



where «Si2(r) = 3<?i • f a% ■ r — a\ • &2- Asymptotically, the radial functions take the 
form 



lim uAr) 

r—*oo 

lim Wd(r) 



-yr 



fl + — 3 

^ d \ 7r (77") 5 



(25) 



where 7 is the wave-number, 7 = ^rriNBd, with Bd the binding energy; rjd is the 
deuteron D/S ratio, and the normalization factor As is defined by 



poo 

A% j o dr 



u 2 d (r) + w 2 d {r) 



(26) 



The experimental values (van der Leun & Alderliesten |1982| ; Ericson & Rosa-Clot 
T983t Rodning & Knutson [T990D are 



Bf p = 2.224575(9) MeV, = 0.8846(8) fm^, rf d xp = 0.0256(4). (27) 
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It is also interesting to note that As can be expressed in terms of the effective range 
defined at the pole-position, 



Pd 



2 / dr 
Jo 



or 



A 2 



u 2 d(r) + wj(r) 
2 7 



With eq. (|28|) , the experimental data eq.(^7]) lead to 



p fP = 1.7635(46) fm. 



7 (1 + V 2 d )A 2 



(28) 
(29) 

(30) 



One can get an idea as to how accurately the Argonne wavefunction describes 
the properties of the deuteron by looking at the predictions: B d = 2.22460 MeV, 
A s = 0.88506 fm^, rj d = 0.02504 and p v d 18 = 1.7660 fm. The close agreement 
between the Argonne results and the experimental values invites us to use the 
two interchangeably in referring to "experiments". Thus for instance in (Park et 
al. |1997b|) , we called the single-particle Ml matrix element given by the Argonne 
wavefunction as "experimental" although strictly speaking it cannot be so called. In 
what follows, we will use eqs.(|2"T|) and ( |3"0"| ) in calculating the Gamow- Teller matrix 
element with the Argonne wavefunction. 



4.2. The pp 1 Sq wavefunction 



An accurate description of the pp wavefunction involves certain subtlety and is 
interesting on its own right. KB (|1994|) made an extensive analysis of the pp wave 
function. Since our result is a bit different from theirs, we discuss here our approach 
to this problem in some detail. 



4-2.1. Potential 



We begin by decomposing the full potential Vf u u into four different components: 
the nuclear potential Vjv, the pure Coulomb potential Vc, the two-photon-exchange 
Coulomb potential Vc2, and the vacuum-polarization potential Vyp: 



= V N + V C + V C2 + V VP 



(31) 



13 



where (Carlson et al. 1991|) 



V c (r) 



a: 



a: 



a 



1 + 2p 2 /m 2 



1 + p 2 /m 2 



(32) 



V C 2 (r) 



Vvp(r) 



(V 2 +p 2 )- + -(V 2 + p 2 ) 



a 



m p r 
2aa' 



Vi 



full\ 



3irr Ji 



dxe 



-2m e rx 



1 A Vx 2 - 1 



2x' 



x- 



(33) 
(34) 



where, in arriving at the second line of eq.([£J), we have removed the non-local piece 
(V 2 + p 2 ) using the equation of motion 



(V 2 +p 2 )4>p P (f) = m p V fu u{r)i) pp (r). 



(35) 



The nuclear potential contains not only strong interactions but also other 
electromagnetic interactions such that the Argonne potential for the pp channel is 
of the same form as eq . (|3l|) .p^| In what follows, when we refer to the "Argonne 
potential" for this channel, we mean eq. fl3~T|) . 



4.2.2. Solution 



Our treatment of the potential V > = Vc2 + Vyp, which is somewhat different 



from that in KB ( 1994 ), goes as follows. First we solve exactly the Schrodinger 



equation with the Coulomb plus nuclear potential 

rv 

V c+N (r) = V c (r) + V N (r) = - + V N (r) (36) 

where we have set a' = a since we are confining ourselves to the case p ~ 0. Then 
V>, which is of order of 0(a 2 ), is treated perturbatively. 

Focusing on the S'-wave, let «c+n be a solution of the radial Schrodinger 
equation for the potential eq.(|36D satisfying the boundary condition m c+n (0) = 0. 



#9 Vc2 an d Vvp i n (Wiringa et al. 1995| ) contain the "form factor" F c (r). We leave out this form 



factor since in our treatment we encounter no singularity associated with these potentials. This 
difference causes little difference in numerical results. 



-14- 



Thus 



dr 2 



+ p 2 - m p Vb+N(V) 



M C +N(r) = 0. 



Let n ~ (10 ~ 20) fm be the range of the nuclear potential so that 

VzvO") =0, r > ri. 



(37) 



(38) 



For r > 7*1, the solution «c+n is given by a linear combination of the regular and 
irregular Coulomb wavef unctions, F and G : 



uc+u(r) = <f>c{r), r>7*i, 

0c(r) = C (77) fG (r)+cot5 c F (r) 



(39) 



with 60(77) = e 2^-\ i V = ^r- Equation ( |39| ) defines the normalization of the 
wavefunction and the phase shift 8 . The properties of the Coulomb functions imply 



0c(O) 
lim 4>c{r) 



Co(v) . 
- — ^ sm 
sm o° 



pr — 7] ln(2pr) + a c + 5 



c 



(40) 



where oq = Argr(l + 777) is the Coulomb phase. The scattering length a c is defined 
by 

\imC 2 ( V )pcot5 c = -^, (41) 
with which the asymptotic wavefunction at threshold reads 



lim0c(r) 



- ^M r ) 



(42) 



with 



0i (r J 
2 (r) 



limC (7 / )G (r) = 2^K 1 (2^), 



(43) 



where i? = (m p a) _1 = 28.8198 fm is the so-called proton Bohr radius, and I v {z) and 
K u (z) are the modified Bessel functions of order v. The scattering length a c can be 
determined by matching the logarithmic derivative at r = 7*1, 



R_ 

Tic 



u 'c+N( r )M r ) - uc+n(r)fl 2 (r) 



and the effective range r c is given by the standard effective range formula 



dr \0" c {r) - u" c+N (r 



(44) 



(45) 
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For the Argonne nuclear potential, the scattering length and the effective range 
come out to be 

a c = -7.8202 fm, r c = 2.782 fm (46) 
in very good agreement with the "experimental values" 

<4 P = -7.8196(26) fm, rf xp = 2.790(14) fm, (47) 



given by the Nijmegen multi-energy analysis (Bergervoet et al. |1988 



In taking into account the effects of the residual long-range potential, 
V> = Vc2 + Vyp, there appear some subtleties. First, although the magnitudes of 
Vc2 and Vyp are small [~ 0(a 2 )}, they are both extremely long-ranged, making 
it difficult to solve the Schrodinger equation directly. We bypass this difficulty by 
treating them perturbatively. Second, depending on how the potential V> is treated, 
there exist in the literature (Bergervoet et al. pL988|) a number of different definitions 
of phase shifts and low-energy scattering parameters extracted therefrom. The 5 C in 
eq.(|39"D represents the phase shift due to Vq+n relative to the phase corresponding 
to the pure Coulomb potential , p iU | Third, we need to keep track of the changes in 
the large-r behavior and in the wavefunction normalization factor which may cause 
difference in the flux density. They should be taken into account self-consistently 



in the definition of the cross section, as pointed out by KB ( |1994|) . We resolve this 
problem in such a way that the V> potential affects the cross section only through 
the reduced matrix element, i.e., by requiring that u pp (r) - the radial wavefunction 
from the full potential - have, modulo change in the phase shift, exactly the same 
large-r behavior as that of «c+n(^): 



hm Upp[r) = — — — sin 



sin 



pr -7] m(2pr) + a c + 5 C + A (4 



where Aq is the improved Coulomb- Foldy correction (Bergervoet et al. |1988j ). This 



means that in order to have the correct flux density, eq. (|48|) requires the full pp 



* 10 One may want to use the phase shift due to the full potential relative to the phase of the 
Coulomb potential, but such a phase shift is inadequate for our present purpose because it results 
in an absurdly big scattering length due to the extremely long range of V>. On the other hand, 
it is possible to use the phase shift due to the full potential relative to the phase corresponding to 
the Coulomb plus V> potential. We shall not adopt this definition because it involves a much more 
complicated effective range formula. 
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wavefunction to be of the form, 



N - 3 S; -° ~ CMaC - (49) 



The u pp (r) that satisfies the boundary condition eq. ([48l) is 

— sin 
u pp (r) = cosAoMc+N(r) 2 [u c +N{r)P v (r) +v c+ N(r)P u (r)] 



P v (r) = m p J dr / fc+N(r / )V>(r / )u P p(r'), 

P u (r) = m p [ <ir / uc+N(r / )V>(r / )Mpp(r') (50) 



o 



where t>c+N( r ) is a second solution of eq.([37]) that satisfies the boundary condition 



vc+x(r) = C (rj) [cos 5 c G (r) - sin5 c F (r)] , r > r\. (51) 

The first order correction in V> can be obtained by replacing u pp by mc+n in P u 
and P v . The threshold limit is readily obtained by noticing that A — > and 
' U c+N( r ; r > ri) — ► 0i (r) when p — >■ 0. 

A remark is in order here about our treatment of the contribution from the C2 
potential given in eq.(^). Up to first order in Vc2 and Vyp, we have 

Vb(r) « (- + VJv(r)) = V c c 2 (r) + V^(r). (52) 

Since the first term, Vq 2 (t), is attractive, its contribution to the phase shift adds 
destructively to the effect of the repulsive VP potential. An approximation frequently 
used in calculating the phase shifts is to take only this term, ignoring the second 
term, Vq 2 (t). This is reasonable for scattering amplitudes since, as one can see from 
eq.flljDD, V c ^(r) cannot affect significantly the large-r behavior of u pp {r). For the 
calculation in question, however, we need to keep the full potential eq. fl52"D . The 
reason is that we are concerned here with not only the phase shifts but also the 
wavefunction in short and intermediate ranges. In fact, we find that the effect of 
Vq 2 (t) is opposite in sign to, and greater than, that of Vq 2 {t) for r <; 49 fm (except 
for very short distances, say, r <^ 0.05 fm). Consequently, the C2 contribution to 
the pp amplitude adds constructively to the VP contribution. These features are 
exhibited in Fig. 1, where we plot the ratio {u pp — uc+n)/mc+n as a function of r. 
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5. Numerical Results 

It is convenient to define the reduced matrix element of the weak current as 

M=M ib + M 2 b (53) 



such that 

(ipd\J~\ip PP ) (~ . A*v 



8irg A A s N pp \ 2g A m N 



k x s d \ M 1B + s d M 2 B 



= iys d -i^^kxs d ^ M + 0{Q\Q 2 \kUk\ 2 ), (54) 

where s d is the spin polarization vector of the deuteron and the normalization factors 
As and N pp have been defined in eqs.(p6|) and (f|jj). Although the reshuffling of 
terms in the last line introduces a term of 0(Q 3 \k\) (i.e., the product of the WM 
term and two-body term), the error made thereby should be negligible. We give our 
numerical results in the increasing orders in chiral counting. 



5.1. Oil) &z, 0{Q) contributions: The single-particle matrix element 



It follows from eqs.([T^) and (0) that the single-particle (reduced) matrix 
element simplifies to 

/•OO 

Mm = / dru d (r)u pp (r). (55) 
Jo 

Following the discussion given in the preceding section, we formally separate the 
contributions into the "Coulomb plus nuclear part" and the "large-r part" 

A^ib = M^ N + M> B , (56) 

/•oo 

M^ N = / dru d (r)u c+ N(r), (57) 
jo 

M m = / dru d (r) [Upp(r) - u c+ ^(r)] . (58) 



Clearly the most important term is the "Coulomb plus nuclear term" eq.( 
so it needs to be calculated as accurately as possible. This term can in turn be 
decomposed into a term which can be extracted from experiments using the effective 
range formula and the remainder that requires theoretical input from the Argonne 
potential: 

MC+N = M ER + ^C+N (5Q) 
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If the next-to-leading-order approximation is made to the effective range formula, 
the resulting matrix element is 



;W 1B 



dr i 



>d{r)<Pc{r) 



+ 



2 Jo 



dr 



1 + 



u 2 d (r) +wj(r) 



R 



where ( = (jR) 1 and 



MO 



1 + ^ 

Ei(C) - 



OO g" 

dt 



- 0Ar 



+ 



«C+nM ~ 4>U r ) 



e^' 
-C* 



_l_ r 



c 



t 



oo p- 

dt 



C 



(60) 
(61) 

(62) 



Now we can evaluate with the experimental values for the low-energy constants 
PD and (0). The result is 



given m eqs. 



cxp 



5.986(1) fm- 1.138(5) fm 



(63) 



where we have given the individual contributions of the two terms in eq.(^T|). If one 
calculates with the Argonne potential instead of using the experimental data, one 
obtains 



M 



ER 
IB 



5.986 fm - 1.137 fm 



Argonne 

which shows once more that it is reasonable to identify the results given by the 
Argonne potential with "experiments". We can now estimate accurately the 
remaining small contribution coming from the difference between eqs.fl57|) and (|60| 
using the Argonne wavefunctions: 



(64) 



5M^ 



0.011 fm. 



(65) 



Since this quantity is very small and the Argonne wavefunction is accurate enough, 
we have not included the error associated with this correction. Combining eqs.( 
and (|65|), we obtain 



IB 



(1 =f 0.02 % =F 0.07 % T 0.02 %) ■ 4.859 fm, 



(66) 



where the first error comes from the uncertainty in a c , the second from r c and the 
third from p d . (The use of "=f" instead of "±" indicates that A^f^ N decreases as 
each of \a c \, r c and pd increases.) 

Finally the effects of the residual C2 and VP potentials are calculated from 
eq.©: 



IB 



mX£ + m 



IB 



-0.031 fm 



(67) 
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with M-Yb = — 0.022 fm and M.^ = —0.009 fm. For later purpose, we evaluate the 
ratio 

5m = = C + ^ = -0-63 %, (68) 



1 1B 



with 



(JVP = _o.45%, 

= + «$<j* N = (0.03 %) + (-0.21 %) = -0.18 %, ((39) 



where we have presented each contribution from Vq 2 and V^. It is worth noting 
that the contribution from the C2 potential is about 40 % of the VP contribution, 
and has the same sign. This is contrary to what one would expect if the Vq 2 were 
ignored. 



We follow KB ( p.994p and express the result obtained up to this point in terms 
of A defined by 

A _ |V^vy AsMm (70) 

In terms of A, our result can be summarized as [see eq.flBDDl p 11 ! 

A 2 = (1 + 5 lB ) 2 (1 ± 0.01 % =F 0.07% ± 0.07%) 2 x 7.02, (71) 
= (1 ±0.01% =F 0.07% ±0.07 %) 2 x 6.93 (72) 

where again the errors correspond to the uncertainties in \a c \, r c and pd (or As) 
in this order. The only possible sources of theoretical error, which are due to 
uncertainties in the Argonne potential, are in dM.^^ (eq.fpSf) and 5pg :N (eq.[p5[]). 
If we were to assign conservatively a ~ 50% error to the two quantities there, the 
resulting theoretical error on the reduced matrix element would then be ~ 0.1 %. 



5.2. Calculation of in effective field theory 



In Section |2|, we argued that the calculation of the single-particle matrix 
element of the electroweak current using the Argonne wavefunctions is equivalent to 



# n If one calculates with the Argonne wavefunctions instead of using the experimental values, one 

finds A 2 \ A rgonne = (1 + ^1b) 2 X 7.03. 
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a "first-principle" calculation using effective field theory of QCD. How this comes out 
in the case of the np capture was discussed in (Park et al. |1997b|) . Here we briefly 
report the result of a first-principle calculation for the leading Gamow- Teller matrix 
element eq.(^). Such a calculation involving the Coulomb potential is interesting 
on its own merit, so we shall describe it in detail elsewhere (Park et al. 1998|) . 



As in (Park et al. |1997b| ), we shall integrate out all fields other than the matter 
field for the nucleon and treat to next-to-leading order (NLO). The potential has 
then the contact interactions used for the np channel plus the Coulomb interaction 



V(q) 



An 
m N 
3 



Co + (C 2 8 ij + D 2< 7 iJ >Vl + Z±Z 2 ^, 



a 



<A<A 



5 ij _ _ 



(73) 



In our case, of course, Z\Z 2 = 1 for the pp channel and for the np channel. Since 
the Coulomb interaction is not separable, we use the Gaussian regulator as 



V{r) 



(2ny 



S A (q 2 )e^V(q) 



with 



exp 



2A 2 J 



(74) 



(75) 



where A is the cut-off. The coefficients Cq )2 for the spin singlet and triplet channels 
and D 2 are fixed by np scattering and the properties of the deuteron. The cut-off A 
defines a regime in which the effective theory is applicable. It is not a free parameter. 
By the rule of (cut-off) effective field theory, physical quantities should not depend 
sensitively on the exact form of the cut-off nor on the exact value of A. Since the 
pion is integrated out, the cut-off should be ~ m n . Indeed as we found in (Park 



et al. |1997b|) , all the static properties of the deuteron are well reproduced for the 
range of cut-off, (140 ~ 200) MeV. p^ For this range of cut-off, the one-body matrix 
element for the np capture comes out to be 4.00 ~ 3.96 which agrees well with the 
result of the Argonne wavefunction ("experiment"), 3.98. The same formalism turns 
out to give 



A^S" N (EFT) = 4.89 ~ 4.85. 



(76) 



# 12 The numerical values of A figuring here could very well be slightly different from the value in 



(Park et al. f 997b ) because of the slightly different regularizations involved. The range of A found 
there to give correct static properties of the deuteron was about (160 ~ 220) MeV. 
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As in the case of the np capture, this agrees closely with what we might refer to 
as the "experimental" value, 4.86, with an error of about 0.5 %. This validates 
our assertion made in Section 3- F i3 | In our approach, the corrections from the V> 



potential are found to be 

«JVP = -0.45 %, 5?£ c = 0.03 %, <fg sN = -(0.19 ~ 0.15) %. (77) 



These values validate the former calculation, eq. (|5"9"|) , where we used Argonne 
wavefunction corresponding to the all order of contributions. Both and S^ :C 
are found to be extremely insensitive to the cut-off. Even £pB :N turns out to be only 
mildly dependent on A. This is somewhat surprising, for the effect of Vq 2 oc ^V]v(r) 
is expected to be amplified at small r by the nuclear potential. 



5.3. 0(Q 3 ) contribution: The exchange-current matrix element 



The matrix element of the two-body current eq.(^) can be separated into 
three pieces, 

M 2B = M^ tc + + M 5 2B (76 

where M 2B ls the contribution from the "kinetic term", Ai 2B i s the zero-ranged 
contributions, and A^b* 6 is the remaining finite-range contribution. From eqs. 
and (H) 



M 



finite 
2B 



A/f kin 
Jvl 2B 



M 5 2B 



2 / 1\ w? r°° 

- [c 3 + 2c 4 + -J ^Jj 2 dru d (r)y Q (m n r)u pp (r) 



2\/2 / A 1 \ III' 



dr w d (r)y 2 (m n r)u pp (r), 



°° dr yi ^ mnr 



3m N f% Jo 2m\r 
y/2ml f°° yi(m n r) 



4/ m N fl Jo 

u<i(r)u'(r) - u' d (r)upp(r) 



(79) 



3mjv/^ Jo 2m\r 
-2 di- 2d 2 + -c 3 + -c 4 + 



w d (r)u' pp (r) - w' d {r)u pp {r) - -w d {r)u pp {r) 



1\ 1 



(80) 



3 3 6/ m N fl Jo 



oo 



S(r) 



dru d (r)^—^u pp (r). (81) 



# 13 While small in magnitude from the view point of the next-to-leading order EFT, the error bar 
of about 0.5 % is bigger than that with the Argonne potential, ~ 0.1 %. This difference may be 
accounted for by the next-order contribution, namely, the effective volume term, which we have not 
considered. 
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It is straightforward to calculate the matrix elements with the Argonne 
wavefunctions. The results are 

M^ te = 0.246 fm, (82) 
M 2B = 0.007 fm, (83) 
M S 2B = -(Ji - 2d 2 + 0.78) x (0.001 fm). (84) 

There are two important points to remark here. The first is that the finite-range 
term eq.(j8~2"|) receives a dominant contribution (~ 93%) from the deuteron D state. 
This means that the matrix element is somewhat sensitive to the shorter-range 
part of the current, i.e., Di^m^r). We will come back to this later. The second 
point is that the zero-range contribution eq.(|84]) is strongly suppressed insofar as 
the magnitudes of the (unknown) constants are "natural", that is 0(1), as 
required by chiral symmetry. Indeed in the way formulated in (Park et al. |1995| ), 
the zero-range operators should be killed by the short-range correlation associated 
with the very short- distance physics (i.e., degree of freedom) that is "integrated 
out" from the low-energy effective chiral Lagrangian. We will therefore drop the 



zero-range term eq.(|84]) and evaluate the finite range terms eq.(|8~2"D and eq.(|33"l 



with a short-range correlation function incorporated as in (Park et al. |1995|) . As 



a short-range correlation function we use a radial wavefunction cut-off of the form 
9{r — r c ) with r c ~ (2m 7r )~ 1 ~ 0.7 fm. 



In Fig. |] is shown the ratio 



* 2B = (85) 



as a function of r c . Taking 0.55 fm < rc < 0.8 fm, our results may be summarized 

as 

5 2B = (4.0 ± 0.5) %. (86) 

This is in the range of the values quoted in the literature (Blin-Stoyle & Papageorgiou 
rgBBJ Gari & Huffman [T9T2] ; Dautry, Rho, & Riska |T9T61 ; Gari |TW§). We also 



present the results with other types of short-range correlation function g[r) that are 



frequently used in the literature (Towner 1987 ) : 



5 2B = 4.84 %, g(r) = 1 - 3o (q c r) = 1 - ( 87 ) 

q c r 

5 2B = 4.01%, g(r) = l-e- ar \l-br 2 ), (88) 
where q c = 3.93 fm" 1 , a = 1.1 fm -2 and b = 0.68 fm -2 . 
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As mentioned, the non-negligible ro dependence in the region rc ^ (2m,,-) -1 
comes from the fact that a short-range (or high-order) current and the D-state 
component of the deuteron wavefunction are involved. This is due to the fact 
that because of the symmetry involved, 0(Q) and 0(Q 2 ) terms relative to the 
single-particle operator are suppressed, an aspect closely related to the chiral 
filter mechanism. This means that some terms higher order than the 0(Q 3 ) ones 
calculated here - starting at 0(Q 4 ) and probing shorter length scales - may not be 
negligible compared with the terms that are retained 1 * 14 



5.4. Cross section 



For comparison with KB (|1994|) , we focus on the low-energy cross-section 
factor Spp(E) (Bahcall and May |1969| ), 

Sp P (E) = a(E)Ee 2nv (89) 

where cr(E) is the cross section for the process ([!]) and rj = Explicitly, it is of 

the form 2 

Spp(E) = (1 + 5 2B ) 2 - m p a G 2 v g 2 A m 5 J(E + E), (90) 
71 7 



where f(Eo + E) with Eq = 2m p — irid is the phase volume defined in (Bahcall |1966 ) 



with however the WM term contribution as well as the deuteron recoil taken into 
account. These two effects turn out to be very small.p^ Setting E ~ 0, we have 



S* PT (0) = 4.05 ( i±*EV ( (*?-)x (10- 25 MeV-barn) (91) 
pp K J \ 1.04 J V1.2601/ V6-93/ V ) \ ) 

where we have substituted Gy = 1.1356 x 10 -5 GeV~ 2 and the phase volume 
f(Eo) = 0.1421. It is worth noting, in passing, that our result for the above S"-factor 



# 14 For instance, if one adopts a phenomenological approach and calculates /9-meson-exchange 
diagrams involving a pNA vertex (with form factors appended to all the vertices), the one-pion- 
exchange contribution is found to be considerably quenched (Bargholtz 1979| ; Carlson et al. 1991). 



From the chiral perturbation point of view, however, these diagrams come at 0(Q 4 ) or higher, and 
therefore they need to be calculated together with all other diagrams of the same chiral order. We 
will not attempt such a calculation here. 

# 15 The WM term increases / by 0.02 % while the deuteron recoil correction decreases it by 0.06 %, 
so the net effect is a 0.04 % decrease. 
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is a precise prediction based upon the Standard Model, that is, without neutrino 



mass (oscillation). The corresponding expression given by KB (|1994|) is 



0» = 3.89 (i±*)' (^) x (10- MeV-baxn). (92) 

By plugging in the numbers obtained above, we arrive at 

S*f T (0) = 4.05 x KT 25 (1 ± 0.15 % ± 0.48 %f MeV-barn (93) 

where the first error arises from the one-body matrix element and the second from 
the two-body term due to the uncertainty in the short-range correlation. (The latter 
does not include the possible contribution of the next order (i.e., 0(Q 4 )) term; see 
later.) This is to be compared with that of KB 

5™(0) = 3.89 x 10~ 25 (1 ± 1.1 %) MeV-barn. (94) 

We are now ready to make a detailed comparison between our result and that 
of KB: 

1. The weak-interaction constants used by KB are: qa = 1.2575 and 

G v = 1.151 x 10" 5 GeV" 2 (Freedman ggg^B Ours are g A = 1.2601 and 
G v = 1.136 x 10~ 5 GeV~ 2 . 

2. KB gave A 2 = 6.92 and §23 = 1 % while we have obtained A 2 = 6.93 and 
S 2B = 4 %. 

Apart from small differences in the constants used, the major difference is in the 
exchange- current contribution. In the next section, we will discuss whether this is 
a genuine difference and to what extent our calculation will be affected by higher 
chiral-order terms. 

Although, as noted by KB ( |1994j ), the corrections due to the long-range 
potential V > are small, it may be of some interest to make a few remarks on its 
role. Our result is that the VP potential decreases the pp rate by about 0.9 % in 
agreement with (Bohannon & Heller |1977|) (where the decrease was found to be 
about 0.8 ~ 1.2 %) as well as with KB (119941) , who found a 0.9 % decrease^. Our 



# 16 The value of Gv is not given explicitly in KB ( |l994|) , but we can reproduce eq.(|9~i|) only with 
the use of this value of Gv , a value which was also adopted in (Carlson |1991| ) . 

# 17 As mentioned, despite the similarity in the numerical results, our calculational framework is 



different from that of KB. This difference will be analyzed in detail elsewhere (Park et al. 1998) 
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result is slightly at variance with that of (Gould |1990| ), wherein with the use of 
the WKB approximation a 1.3 % decrease was found. In addition the C2 potential 
lowers the pp rate by about 0.4 %. In total, therefore, the downward correction due 
to V> is somewhat larger than that of KB. Since the A 2 is about the same, this 
means that our single-particle matrix element Aim is somewhat larger than that of 
KB. 



6. Discussion 

Following the procedure of chiral perturbation theory proven to be highly 
successful for the thermal np capture process, we have calculated the pp fusion rate 
to 0(Q 3 ) in chiral counting relative to the leading single-particle Gamow- Teller 
matrix element. To the order considered, the error involved in the calculation is 
small, <^ 1 %. This result, given a justification from a cut-off effective field theory of 
low-energy QCD as in the case of the np capture (Park et al. |1997b|) , supports the 



canonical value of Bahcall et al. and does not support the "relativistic field theory 
model" result of Ivanov et al. (|1997| ). 

The main caveat in this calculation is in the meson-exchange contribution 
which comes out to be about 4 % when calculated to the chiral order 0(Q 3 ). At 
the next order, 0(Q 4 ), loops and higher-order counter terms enter, so that there 
is no reason to believe that they are negligible compared with the 0(Q 3 ) tree 
contributions. (For instance, it could be lowered to about 1 ~ 2 % instead of ~ 4 % 
found here). This aspect is different from the case of the np capture where the chiral 
filter mechanism assures the dominance of the tree-order pion exchange-current 
contribution. Here the absence of the chiral filter phenomenon can allow higher-order 
(loop) terms to figure equally importantly as the tree-order terms. From this 
viewpoint it is not surprising that a model calculation of the terms of 0(Q 4 ) and 
higher based on the vector-meson exchange and form factors (Bargholtz |1979| ) 
indicates that there can be a considerable suppression of the tree-order correction. 
Although such a reduction in the exchange- current contribution seems to go in the 
right direction for beta decays of higher- mass nuclei (Carlson |1991|) , it is probably 



unsafe to import the result of such model calculations into our work. To go to 0(Q 4 ) 
or above in our theory at which heavy mesons and form factors come in, a large 
number of Feynman graphs of the same chiral order have to be computed on the 
same footing to assure chiral symmetry, and such a calculation has not been done 
yet. In the absence of consistent calculations, our attitude is that we should not 
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attach any error estimates on terms not accounted for to the chiral order computed. 
Calculating the higher-order terms will be left for a future exercise. 
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r (fm) 

Fig. 1. — The distortion of the pp wavefunction due to large-r EM residual interactions 
defined as (u pp — uc+n)/uc+n is shown for various components entering into u pp . The solid 
line denotes the distortion due to the total (VP+C2) residual interactions, the dotted line 
due to the VP, and the dashed line due to the C2. The contribution from C2 is decomposed 
into that due to V§ 2 (dot-dashed) and that due to Vq 2 (dot-dot-dashed), see eq.(|52]). The 
C2 contribution changes from negative to positive at r ~ 49 fm. 
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